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Abstract 



We study the decay B — > K£ + £ for £ = e, u,t with a softly recoiling kaon, that is, for high 

i 

O , dilepton invariant masses \fq* of the order of the 6-quark mass. This kinematic region can be 

<D 

treated within an operator product expansion and simplified using heavy quark symmetry, leading 

^ to systematic predictions for heavy-to-light processes such as B — > R(*H + £~ . We show that the 

> ______ 

decay rates of both B — > K*£ + £ and B — > K£ + £ decays into light leptons depend on a common 

combination of short-distance coefficients. The corresponding CP-asymmetries are hence identical. 

(N 

Furthermore we present low recoil predictions for B — > K£ + £ observables, including the flat term 
-_H in the angular distribution which becomes sizable for taus. We work out model-independently the 

constraints on AB = 1 operators using the most recent data from the experiments BaBar, Belle, 

•i-H 

^ CDF and LBCb. For constructive interference with the standard model, generic new physics is 

pushed up to scales above 44 TeV at 95% CL. Assuming none or small CP-violation we obtain a 
lower bound on the position of the zero of the forward-backward asymmetry of -B — > K*°£ + £~ 
decays as > 1.7 GeV 2 , which improves to q$ > 2.6 GeV 2 for a standard model-like sign b — > 57 
amplitude. 
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I. INTRODUCTION 



Rare 5-decays into dileptons are precision probes of the standard model (SM) and the fla- 
vor sector and provide constraints on physics beyond the standard model (BSM). Important 
semileptonic modes in terms of experimental accessibility and theory control are those into 
a if or a K* . The latter decays exhibit a rich phenomenology especially through angular 
analysis of subsequent decays B — > K*{-^- Ktt)£ + £~, see [1] for a recent summary. Decays 
into a pseudo-scalar meson B — > K£ + £~, £ = e,fi,T allow to perform a number of comple- 
mentary measurements as well. Observables include the decay rate 1^, its CP-asymmetry 
A l CP , the forward-backward asymmetry A FB and the flat term Ffj appearing in the angular 
distribution [2HD 

j 3 (l-F£)(l-cos 2 ^) + ^ + 4 B cos^, (1.1) 



or distributions in the dilepton mass thereof. (The angle 9^ is defined in Section II C ) 
Previous systematic analyses of B — > K£ + £~ distributions [4J focused on the region of 
large hadronic recoil, where QCD factorization (QCDF) applies [51 E]- The intermediate 
recoil region where charmonium-resonances dominate the dilepton spectrum through B — > 
K(J/^f,^>') — > K£ + £~ decays has been studied recently [7j. In this work we provide a 
systematic analysis for the region of low hadronic recoil, that is, for large dilepton masses 
\f(f of the order of the 6-quark mass, above the \l/'-peak. 

The study of heavy-to-light decays at low recoil by means of a local operator product expan- 
sion (OPE) in has been put forward by [S], and recently see also [TU] for earlier 
mention on inclusive decays. The simultaneous matching onto heavy quark effective theory 
(HQET) and use of the improved Isgur-Wise form factor relations [11] has been shown to 
be of benefit for B — > K*£ + £~ decays by allowing for the design of specific low recoil observ- 
ables with sensitivity to either short- or long-distance physics as well as checks of the theory 
framework \1'2\ [T3]. Equally important, the theoretical accessibility of this kinematic region 
is necessary for a full exploitation of the available and future rare 5-decay data. Note that 
roughly a similar amount of B — > K*£ + £~ data exists presently for the low and the large 
recoil region by all contemporary 5-physics experiments, BaBar, Belle, CDF and LHCb. 
The latter will soon be competitive and eventually take over in statistics as indicated by the 
recent preliminary results on the angular distribution of B° — > K*°£ + £~ decays [T¥l 15J. As 



for B — > K£ + £ decays, LHCb reported so far 35 ± 7 events of B + — > K + fi + fi within a 
recorded luminosity of 37pb _1 during the 2010 LHC run |16j . 

The plan of this paper is to work out and explore the phenomenology of B — > K£ + £~ decays 
at low hadronic recoil. We give in Section [TT] the heavy-quark form factor relation, collect 
the expressions for the decay amplitude at low recoil and discuss the observables relevant 



to our model-independent framework. In Section III A we present SM predictions including 
a discussion of their uncertainties. In Section IIII Bl we derive constraints on the AB = 1 
short-distance couplings entering B — > K^£ + £~ decays using the most recent experimental 
low and large recoil data. We conclude in Section |IV| Details on the form factors are given 
in Appendix [A} 



II. THE DECAY B -> Ki + l~ AT LOW RECOIL 



We use an effective AB = 1 Hamiltonian to describe the flavor-changing b — > s£ + £ transitions 

as 

%* = -%-v»K (2.1) 

a/2 k 

x (C 7 ^ [sa, u P R b] + C 9 [s^P L b] [hH] + C w [s^P L b] [l^d] + h.c.) + . . . . 

Here, the ellipses denote tree-level induced or subdominant contributions which we assume 
to be SM-like. Furthermore, Vij denote CKM-elements and nib the MS-mass of the 6-quark. 
We allow the Wilson coefficients C 7 9 10 to be complex- valued to account for CP-violation 
beyond the SM. We assume them throughout this work to be evaluated at the scale /x = m b . 
For further details we refer to previous low recoil works employing the same notation [T21 [T5] . 



A. The Improved Isgur-Wise Relation 

The matrix elements of B — > K transitions can be parameterized in terms of three q 2 - 
dependent form factors /+,t,o which are defined in Appendix |Aj Following [HI [TTJ , the QCD 
operator identity 

id v {si(j^b) = id^{sb) - m b (s7 M 6) - 2(siD^b) (2.2) 
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allows to derive an improved Isgur-Wise relation between fx and /. 



+ • 
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(2.3) 



(2.4) 



in agreement with [17] . (We neglect the mass of the strange quark.) Here, we denote 
by m B) rriK the 5-meson and kaon mass, respectively, and neglect in the second line the 
subleading HQET form factor 5+\ see Appendix [X The l/g 2 -factor on the right-hand side 
of Eqs. (2.3) and (2.4) is of kinematical origin, related to the definition of fx, Eq. (Al ). The 
/x-dependent coefficient k reads, including 0(a s ) corrections, as 

, D o\v)\ ™> b {fj) 



m B 



with the HQET Wilson coefficients C^\d^ given in [8j. At /i — m& holds k — 1 



While deriving Eq. (2.4) one finds that the form factor f is a power correction 

/ A \ 



/ (g 2 ) = o + c)! 

as expected from heavy quark symmetry [18J. 



\m b ) 



(2.5) 
0{al). 

(2.6) 



B. 



The B -»• Kt 



Matrix Element 



Beyond the contributions from the operators in Eq. (2.1 ), the amplitude of B — >■ K£ + £ also 



receives contributions from current-current and QCD-penguin operators. As a result of the 
OPE in 1/ yfq* for the leading dimension 3-operators, these contributions can be taken into 
account by effective Wilson coefficients Cj 9 [9], whereas sub- leading contributions enter at 
dimension 5 and are suppressed by (A/m b ) 2 ~ 2%. Here we follow [8] and subsequent low 
recoil B — > K*£ + £~ works [T2l [TB^ and perform an additional matching onto HQET and use 



the form factor relation Eq. (|2.4|). The effective coefficients are then given as 



Ceff n 
7 — 7 



47T 



Cf = C 9 + h(0,q 2 ) 



h(m b , q 2 



A -c +c 



■C. + 52C, 



32 



(2.7) 
(2.8) 
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7C 3 + -C 4 + 76C 5 + y C 6 
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C x (B(q 2 ) + AC(q 2 )) -3C 2 (2 5(g 2 ) - C(g 2 )) -C 8 F 8 (9) (g 



m 



2 r 



9 Cl + 3 C2 



1 + A U ) + 2C 3 + 20C 5 



These include the NLO QCD matching corrections and doubly Cabibbo-suppressed contri- 
butions proportional to \ u = V ub V* s / (V tb V t * s ). The latter are responsible for the tiny amount 
of CP-violation in the SM in b — > s transitions. 

The B — > K£ + £~ decay amplitude [4] simplifies within the SM operator basis Eq. (2.1) after 



applying the form factor relation Eq. (2.4) to 

.Gpa e 



where 



KP 



V2tt 



F a = C 



10 J 



V tb V* f + (q 2 ) [F vV ^h/) + F A p»(£ W £) + F P (I l5 £)] , (2.9) 
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10 > 



(2.10) 



and denotes the 4-momentum of the I?-meson and mi the lepton mass. 



The use of Eq. (2.4) introduces an uncertainty of 0(A/m b ) in the term proportional to Cj in 
Fy, however, the phenomeno logical impact of this uncertainty is additionally suppressed by 
C^ s /Cg S , which is ~ 0.1 in the SM. Further subleading I/nib contributions to the amplitude 



Fy in Eq. (2.9) itself receive an additional suppression of a s [8]. The associated uncertainties 
are included in our phenomenological analysis following the procedure described in |13j . 
The coefficient Fp is suppressed by mg/iriB and hence small for light leptons. Since the form 
factor /o enters Fp only, its impact for phenomenological implications is negligible except 
for taus. 



C. Observables and Angular Distribution 



Continuing along the lines of [1], we write the differential B — > K£ + £ decay distributions 



as 



d 2 T e [B -+ K£+t 
dq 2 d cos 9 1 

dTf[B ->• K£+t 
dq 2 



? (<? 2 ) + be(l 2 ) cos 9 a + c e (q 2 ) cos 2 9 e 



a i(l 2 ) + 3 C K? 2 ) 



(2.11) 
(2.12) 
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with g 2 -dependent observables ae, hi and q. The angle 9e is defined as the angle between 
the .B-direction and the ^"-direction in the £ + £~ rest frame. Within the SM operator basis 



Eq. (2.1) 



with 



A 



\Fa\ 2 + \F v \ 2 ) 



+ 2m £ (m| - m 2 K + q 2 )Re(F P FX) + Am] m 2 B \F A \ 2 + q 2 \F P \ 2 



0. 



A 



-Plj(\F A \ 2 + \Fv\ 2 ) 



G 2 F a 2 e \V tb V t * s \ 2 
2 9 7r 5 m| ' 



ft 



4mj 
q 2 



(2.13) 



(2.14) 



\ \ / 2 2 2\ 4 i 4,4 rW2 2, 22, 2 2\ 

A = X(m B , m K , q J = m B + m K + q -2 [m B m K + m B q + m K q ) . 



The terms in the B — y K£ + £ angular distribution Eq. (1.1) can then be obtained as [I] 



a 2 

l imax 



dq 2 



1 

— ( 

3 



A FB — 0, 



H 



cig 2 (a£ + a) 



(2.15) 
(2.16) 
(2.17) 



The forward-backward asymmetry of B — > K£ + £~ decays is zero within the SM operator 
basis^j] and Ffj is suppressed by lepton mass, see below. However both observables receive 
contributions from scalar and tensor operators, and can signal such BSM effects [1, 20J. 
For £ = e, /i we can safely neglect the lepton mass. In this limit ft = 1 and 



at = r — — f + pi, 



with the short-distance coefficient 



Pi 



s 2mbmB s 



\c 



in I 



(2.18) 



(2.19) 



1 Non-zero values in the SM originate from QED corrections and are tiny [U [TH] . 
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The differential decay rate at low recoil can then be written as 

dr ' [B : q r n ^ffi^ (2,o) 

and Fjj vanishes for = 0. 

The SM value of p\ and its associated uncertainties have been presented earlier in the 
discussion of B — y K*£ + £~ decays at low recoil [12J. The appearance of the same short- 
distance factor pi in both vector and pseudoscalar final states uniquely correlates the decays 
B — y K£ + £~ and B — y K*£ + £~ into light leptons. As a consequence, the CP-asymmetry of 
the decay rate, A^-p, i s identical 

a" \b > wn-\ - dTi[B ± R ^nW - Wt[B ± Kfe-]W (221) 



Pi- Pi _ (l) 



a^'lB -»■ K 
Pi + Pi ° 

where p~\ is obtained from p\ by complex conjugation of the weak phases, i.e., the CKM 
matrix elements and the Wilson coefficients C i . The CP-asymmetry Oq P is form factor-free 
and essentially vanishes in the SM |a^p|sM < 10~ 4 [T3] . 
Keeping the lepton mass finite we obtain 

To 



n 4m 'ir \ 2 T 
Pi ^2~l L 10l *M> 



Q= 4 



V °V\ 3 p!f + Pu (2.22) 



F = l- {m *- m « )2 ff) . (2.23) 



4 

with 



The second term in Tq is order one despite fo being a power correction, see Eq. (2.6 ), because 
of the phase space enhancement by the factor A ~ (D(A 2 m 2 B ) in the denominator. 
The ^-differential flat term of the angular distribution 

e 2 a e + Cj 6mj Pi ~ l^iol 2 ^" 3 (r> 0A \ 

*h{q ) = _ , i. = — t~ x y z - 2A > 



0) 

2\ 



can be sizable for £ = r, as we show explicitly in Section III A for the SM. The FJj{q 
observable is complementary to the B — y K£ + £~ branching ratios. From combined analysis 
with the branching ratios of B — y K^*H + £~ for £ = e, p decays the coefficient p\ could be 
extracted together with the magnitude of C 10 , see also [21], given sufficient control of fo/ f+- 



A 


n 01 O+0.013 

— U.UZ ( 


123 


A 


0.22543 ± 0.00077 [2H] 


P 


0.144 ±0.025 


M 


V 


0.342+2-SJjj [23] 


oi s {M z ) 


0.11762 




T B+ 


1.638 ps [21] 


a e (mfe) 


1/133 




T go 


1.525 ps El] 


m c (m c ) 


(1.27+Hq) GeV 
v — u.uy/ 


m 

t 3 


mg+ 


5.2792 GeV El] 


rrib(mb) 


(4.19+g-JI) GeV 


M 




5.2795 GeV [S] 


pole 

t 


(173.3 ±1.3) GeV 


[251 

! J 


m K + 


0.494 GeV [2] 


m e 


0.511 MeV 


m 


m K o 


0.498 GeV [M] 


my, 


0.106 GeV 


m 


M w 


(80.399 ± 0.023) GeV [21] 


m T 


1.777 GeV 


m 


sin 2 6 W 0.23116 ± 0.00013 [M] 



TABLE I. The numerical input used in our analysis. We neglect the mass of the strange quark. 
The superscript °( + ) refers to neutral (charged) parameters. 

III. PHENOMENOLOGICAL ANALYSIS 

For the phenomenological analysis of the various rare B decay observables we use the nu- 
merical inputs given in Table [Tj All numerical results presented are obtained with the EOS 
flavor tool [22]. All experimental B — > K^*H + £~ results shown and used are CP-averages, 
and should silently be understood as those. 



A. Standard Model Predictions for B -> K£ + £ 

The SM predictions for the B — > K£ + £~ observables at low recoil in the framework described 
in Section [TT] are given in Table [TT| The integration region is 14.18 GeV 2 < q 2 < (mB—mx) 2 , 
i.e., above the ^/'-resonance. We employ the form factors /+ and fo - the latter is required 
for decays to massive leptons only - from light cone QCD sum rule (LCSR) calculations and 
a g 2 -shape by a series expansion from [7] extrapolated to high q 2 . Our choice is motivated 
by the good agreement with preliminary unquenched lattice results [2S] and the availability 
of uncertainties for the parameterization of the form factors, see Appendix [A] for details. 



The dependence on the tensor form factor has been removed by means of Eq. (2.4). 



The uncertainties given in Table [TT] originate from the form factors (FF), specifically the 



variation of the inputs Table IV to the parameterization Eq. (A3), the CKM matrix elements 



8 



Observable 


Central Value 


FF 


SD 


SL 


CKM 


10 7 x B(B~ -+ K-t+f) 
10 7 x B(B~ -»■ K-T+T-) 

10 3 x 

10 1 x Fj5 


1.04 
1.26 
7.50 
8.90 


+0.60 
—0.27 

+0.40 
-0.21 

+2.85 
-2.61 

+0.33 
-0.45 


+0.04 
— 0.02 

+0.04 
-0.03 

+0.04 
-0.10 

+0.01 
-0.02 


+0.03 
—0.03 

+0.02 
-0.02 

+0.10 
-0.10 

+0.02 
-0.02 


+0.04 
—0.07 

+0.05 
-0.09 



TABLE II. The SM branching ratios and flat terms F B Eq. (2.17) integrated from 14.18 GeV 2 to 
(mg — itik) 2 - For a description of the uncertainties see text. The branching ratios for the light 
leptons £ = e and £ = fj, are equal within the given precision (finite lepton masses are taken into 
account). The branching ratios for neutral B decays are obtained by rescaling the corresponding 
charged ones with t b o /t b + . All branching ratios are identical within the given accuracy to the 
ones of the corresponding CP-conjugated decays. 

(CKM), and the short- distance parameters m^ ole , Mw, sin 2 9w as well as the scale \i varied 
within [/ift/2, 2/ib] with ^ = 4.2 GeV collectively denoted as (SD). We evaluate the SM 
Wilson coefficients at next-to-next-to-leading order. In addition, we take into account the 
(unknown) subleading 1 / m^-corrections to the B — > K£ + £~ matrix elements and form factor 



relation Eq. (2.4) (SL). The latter uncertainty is estimated with the procedure outlined in 
the Appendix of [13] adopted to B — > K£ + £~ decays. The residual CKM uncertainty from 
X u in Ffj has been dropped in Table [n] because it is below the given precision. To the 
accuracy given in Table [TT| SM CP-violation is too small to induce any effect. 
At high q 2 weak annihilation topologies are strongly suppressed by (A/rrib) 3 [9], and thus 
isospin breaking effects are dominated by the differences in the lifetime. The branching 
ratios of decays of neutral I?-mesons are obtained by rescaling the corresponding branching 
ratios of charged B decays in Table |n] with the ratio of lifetimes r B o / r B + ■ In the numerics 
we keep as well the very small effect of isospin breaking by the masses m# and itlk- 
We show the SM predictions for dB(B~ — > K~ /i~) /dq 2 in Fig. [I] versus existing data[^] 
by Belle [27] (red points), CDF [2H] (black points) and BaBar [3] (orange points) and with 
individual uncertainty budgets in Fig. [2j In Fig. [3] we show the differential SM branching 
2 The data on the binned B — > Kt + t~ branching ratios by Belle and BaBar are publicly available for an 
unknown admixture of charged and neutral B decays only. We find differences of the order of the scan 
resolution in the constraints Figs. [5] and [6] when interpreting the Belle data as either purely charged or as 
purely neutral B decays. This effect will be even less important in the future with improved statistics. 
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q 2 [GeV 2 ] 



FIG. 1. The SM branching ratio dB(B~ — } K~ fi + fi~)/dq 2 using form factors from [7] (blue solid 
and dashed lines) alongside the measurements by Belle [27] (red points), CDF [28] (black points) 
and BaBar [3] (orange points), see footnote [2j The blue band shows the theoretical uncertain- 
ties described in the text added in quadrature. The vertical (grey) bands are the experimental 
veto regions [271 28J to remove contributions from B — > J/^(— > fi + [i~)K (left-hand band) and 
B —> > ^i + jjL~)K (right-hand band). The dashed (blue) lines below the J/^f (between the 
charmonium bands) correspond to theory extrapolations from the low recoil (low and large recoil) 
region. 

ratio into ditaus. The evaluation of observables in the large recoil region is carried out within 
the framework of QCDF [51 E] • We employ the same B — > K£ + £~ form factors used at low 
recoil, i.e., the ones from Ref. [7J. The treatment of the subleading corrections for low q 2 is 
adapted from [29]. The flat term F^{q 2 ) for both muons and taus is shown in Fig. |4j The 
growth of F^(q 2 ) at both ends of the spectrum results from a finite muon mass. While at 
low q 2 such effects are kinematically enhanced by m^/q 2 , the effect at the very high g 2 -end 
is induced by the vanishing decay rate in the denominator. For electrons these effects are 
further suppressed by ml/m 2 ~ 10~ 5 , and negligible. 

B. Model-Independent Analysis 

We perform a model-independent analysis in the AB = 1 complex-valued Wilson coefficients 
C k = \C k \ exp(z0fc), k = 7,9, 10. We build on the analysis by [121 H3] using the updated data 
from CDF [281 ED], the first B -± K*fi + fx- data from LHCb HUES], and include B K£ + £~ 
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10 12 14 
q 2 [GeV 2 ] 



FIG. 2. The SM branching ratio dB(B~ — > K~ /i + /x~)/dg 2 with the linearly added uncertainties 
from the form factors (blue), the CKM matrix elements (yellow), the short-distance input (cyan) 
and the subleading l/m& corrections (red), see also Fig. [TJ 
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FIG. 3. The SM branching ratio dB(B — > K t + t )/dq 2 with linearly added uncertainties, see 
Fig. [2] The grey band refers to the cut from BaBar's search for this channel [30j . 



decays. The corresponding experimental results for the latter branching ratios from Belle [27] 
and CDF [2E] are given in Table [TTT| and also shown in Fig. [I] (see also footnote [2]) . The data 
by BaBar [3] shown as well in Fig. [I] are not used as none of the two bins are fully applicable 
to low- or high-g 2 theory frameworks. We discard CDF's data on B(B° — > K°[i + /j,~) [28] 
as the uncertainties are very large. First, binned measurements of Ap B for light leptons are 
available and are consistent with A^ B = [271 13T] . These data are not taken into account 
because within the context of the SM operator basis used in this analysis Ap B vanishes, see 
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FIG. 4. The flat term Ffj(q 2 ) for B — >■ K £ + £ decays given in Eq. (2.24) for I = /i (lower curve) 
and i = t (upper band) in the SM. The theoretical uncertainties are added in quadrature. 



Observable 


SM Prediction 


Measurement 


B(B° -> #V/Oi.o,6.o 
B{B~ -> K-fi+fi-) i.o.6.0 


(1.59^1) x 10- 7 *[4j 
(l.75^)xl0- 7 *[4j 


(1.3618S) x io- 7 ti^ 
(l.41±g;|») x 10- 7 M 


B(fl° -> ^V+ M -) 14 . 18)16 .oo 

B(S _ -)• i^-Ai + /i-)l4.18,16.00 


(OMt^l) x 10- 7 
(0.37^°) x 10- 7 


(0.38i°;^) x IO" 7 t[27| 
(0.53±g;if) x IO" 7 [28] 


£(5° -> K /i + /O 16.00,22.86 
-)• if~/X + //~)l6.00,22.90 


(0.63^;1) x io- 7 
(0.681^) x 10- 7 


(0.98±8;||) x IO" 7 tJ27J 
(0.48lg;J|) x IO -7 [28] 



TABLE III. The B — > Kl + £ branching ratios and their respective measurements with systematic 



and statistical uncertainties added linearly by Belle [27] and CDF [28] used to constrain C 



7,9,10' 



All 



entries are understood as CP-averaged. The integration binning g^ in < q 2 < <?max i n GeV 2 is indi- 
cated as B n 2 _2 • * includes updates to the numerical input and uses a different parametrization 

"mill '"max 

of the form factors than [3]. ^see footnote ^1 



Eq. (2.16). Decays into ditaus have not been observed. The presently best limit stems from 
BaBar B(B + -»■ K+t+t') < 3.3 x 10~ 3 for q 2 > 14.23 GeV 2 at 90% CL [3U]. 

We perform a six- dimensional parameter scan over the magnitudes \C 7g 10 | and weak phases 
07,9,10- The resolution in |C 910 | is 0.25, and tt/16 in 07,9,10- We vary 0.3 < \C 7 \ < 0.4 with 
step-width 0.02. Statistical and systematic experimental uncertainies are added linearly in 
the scan. Contrary to [13] we do not symmetrize the uncertainties, however, we checked 
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\Cg\ vs. |Cio| for all data 
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FIG. 5. Constraints on \C 9 \ and |C 10 | from the combined analysis of B — > {K , K* , X s } £ + £ decays 
at 68% CL (red area) and 95% CL (blue area). The 68% CL (dotted) and 95% CL (solid) contours 
without using B — > K£ + £~ decays are shown as well. The green square marks the SM prediction. 

that with the current data the results exhibit insignificant differences only between using 
symmetrized and un-symmetrized uncertainties. The constraints on |C 9 |-|C 10 | from the joint 
B ->■ K*£ + £-, B -> K£+t and B ->■ XJ + £- analysis are presented in Fig. § The 
inner (outer) colored areas are allowed at 68% CL (95% CL). The inclusion of the decay 
B — > K£ + £~ improves the scan, as can be seen from the overlaid contours obtained without 
using B — > K£ + £~ data. As the latter enter the analysis only via branching ratios, they put 
constraints on pi, roughly ~ |C 9 | 2 + |C 10 | 2 , and show very little sensitivity to CP phases [32J. 
In Fig. [6] we show the constraints from B — > Kl + l~ decays alone in the large recoil region, 
the low recoil region and for both kinematic regions, respectively. All these constraints are 
consistent with each other and the agreement with the SM is good. 
From the full scan we obtain the allowed ranges 



2.3 < \C g \ < 6.5 



;i.o < \c g \ < 7.o) 



(3.1) 
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FIG. 6. The constraints from B — > K£ + £~ decays on |C 9 | and |C 10 | from the large-recoil region 
(left-hand plot), the low-recoil region (plot in the middle) and both (right-hand plot), for further 
notation see Fig. [5] 
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FIG. 7. The constraints from all the currently available data on the real- valued Wilson coefficients 
C 9 and C 10 for SM-like sign C 7 < (left-hand plot) and C 7 > (right-hand plot), for further 
notation see Fig. [5] 



2.3 < |C 10 | < 5.3. 



[1-8 < |C 10 | < 5.5) 



(3.2) 



at 68% CL (95% CL), see Fig. |5] This implies that the branching ratios B(B S -> £+£~) oc 
|C 10 | 2 for £ — e, /Ji,t can be enhanced by at most a factor of 1.7 at 95% CL with respect to 
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FIG. 8. The zero of the forward-backward asymmetry q g in B° — > K*°£ + £~ decays at 1 a versus 
\C g \ for C 7 < (solution A in Fig. [7J red lines) and C 7 > (solution D in Fig.JTJ blue lines) assuming 
no BSM CP-violation. The vertical lines denote the corresponding 95% CL upper limits on |C 9 | 
from Fig. [7| The green area marks the SM (for fixed /i = 4.2 GeV). 

their SM value. 

Furthermore, we obtain constraints for real-valued Wilson coefficients C g 10 by discarding 
from the complex-valued scan all data with ^7,9,10 7^ 0,7r. These constraints are shown in 
Fig. [7j A previously existing region, labelled B, is now excluded, while region C is disfavored 
at 95% CL. With CP- violation neglected the upper bound on \C g \ provides a lower limit on 
the zero crossing q g of the forward-backward asymmetry in B — > K*£ + £~ decays [33]. We 
find <2q > 1.7 GeV 2 , which gets improved to q 2 > 2.6 GeV 2 assuming the sign of C 7 to 
be SM-like, that is, negative. The zero as a function of \C g \ (for fixed /1 = 4.2 GeV) for 
B° — > K*°fi + n~ decays is shown in Fig. [8] for both sign combinations which allow for a zero, 
that is, solution A and D of Fig. [7J There is no such limit on q$ in the fully CP-violating 
case. The dashed curves denote extrapolations to q 2 > 6 GeV 2 where QCDF may not be 
applicable. 

In the SM we find the zero to be located at 



2,SM 



("X Q7+0-03I +0.09 1 +0.29 1 \ p„v 2 
yo.Ji _ .03IfF -0.09lsL -0.27lsD/ ^ V 



(3.3) 



consistent with P,[29j|3l]. The SD uncertainty without the one from the scale /1 is l^os GeV 2 . 



The location of the zero is form factor-independent only at lowest order in 1/mb [21 15]. The 
uncertainty denoted by FF stems from varying the LCSR inputs to the form factors used in 
the scan to describe B — > K*£ + £~ decays by [35J. 
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One can ask about the implications for generic new physics from the AB = 1 BSM operators 
J2i T^Oi, where 



9 = S7 M (1 - 75)6 /U7> , Ova = S7 M (1 - 75)^7 M 75A* • ( 3 - 4 ) 

Assuming unsuppressed contributions, |cg ; io | = 1, the scale of new physics Anp must be as 
high as 

A N p > 30 TeV (15 TeV) from 9 , (3.5) 
A NP > 44 TeV (16 TeV) from O 10 . (3.6) 

On the other hand, with Anp = 1 TeV the coefficient of the higher dimensional operator 
O10 needs a (flavor) suppression as strong as |c 10 1 < 5 x 10~ 4 (4 x 1CT 3 ). The corresponding 
numbers for Og read |cg| < 1 x 10~ 3 (5 x 10~ 3 ). The bounds are obtained at 95% CL from 



Eqs. (3.1) and (3.2). For all cases the first (second) number corresponds to constructive 



(destructive) interference with the SM. The bounds from C 10 are stronger than those of C 9 . 



IV. CONCLUSION 



With event samples of order several hundred analyzed [3j [HI [15], [23, |28] and beyond a 
thousand from LHCb alone at the horizon [16] the studies of rare decays B — > K^£ + £~ 
begin to reach deep into the space of AB = 1 weak operators in terms of the allowed ranges 
and correlations of the short-distance couplings. 

Allowing for BSM CP-violation the current status is shown in Fig. |5j updating the analysis 
of [13] for the moduli of the Wilson coefficients C g and C 10 by including the data from LHCb 
[14"] and B — > K£ + £~ decays. For the latter we worked out predictions for the kinematic 
region of low hadronic recoil, proven already beneficial for B — > K*£ + £~ decays [12]. The low 
recoil framework is equally applicable to other S-meson to pseudo-scalar decays, including 
B n£ + £-, B s K£ + £~ and B s r}^£ + £~. We find that the SM is presently in good 
agreement with the rare decay data. 

Assuming none or negligible CP- violation only the current situation is shown in Fig. [7} The 
preferred confidence region is shared between disconnected islands, one including the SM. 
Only two, solution A and D, allow for a forward-backward asymmetry zero in B° — > K*°£ + £~ 
decays. We predict the location of the zero to be above 1.7 GeV 2 . If a zero is found, the 
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solution C, which is already disfavored, is excluded. Vice versa, if no zero is found, solutions 
A, the SM, and D are excluded. To distinguish the solutions A and D requires that the 
precision of the model-independent analysis reaches the level of being sensitive to interference 
terms with SM dominated four-fermion operators. 

The presence of the pseudo-scalar channel improves the statistical power of the analysis. 
Even stronger advantages arise from correlations among the rare decays: We find that at 
low recoil the B — > K£ + £~ amplitude involves the same short-distance couplings present in 
B — > K*£ + £~ for light leptons (£ = e or fi), and the corresponding CP-asymmetries in the 



rate, a9,p as in Eq. (2.21), are identical. Moreover, the angular distribution in B — > K£ + £ 



is an indicator for operators beyond Eq. (2.1). In this operator basis, Afb(B — > K£ + £ ) 



vanishes and Ffj receives only small corrections of the order rnj/m^, hence vanishes for 
practical purposes for the light leptons as well. These observables will therefore be important 
null tests for effects from BSM operators, complementing the search for B s — > fi + fi~ decays 
and lepton non-universality by comparing B — > K£ + £~ observables into electrons with those 
into muons, e.g., [I]. 



Note added: During the completion of this work a systematic analysis of AB = 1 constraints 
including low- and high-g 2 data onB-f K*£ + £~ but without B — > K£ + £~ decays appeared 
|36j . The findings for the SM operator basis are qualitatively consistent with [T2l 113] and 
the present work. 
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FIG. 9. The B — > K form factors /+, fx and /o by [7] extrapolated to the high-g 2 endpoint and 
preliminary lattice results with statistical uncertainties only (data points) [26]. The shaded bands 
show the respective form factor uncertainties. 

Appendix A: The B — ^ K form factors 

The B — > K form factors /+,t,o are defined as usual 

(K(k)\ rfb \B(p)) = f + (q 2 )(p + ky + [/ (g 2 ) - f + (q 2 )] ^ ~ m * <f , (Al) 
(K{k)\so^b\B{p)) = i \{p + kyq v - q »(p + k) v ] , (A2) 

where p (k) denotes the four-momentum of the 5-meson (kaon) and q = p—k. The numerical 
analysis in this work is performed with the form factors from Ref. [7] parameterized as, 
% = +,T,0 



f i ( s ) = -^^^l + b\iz(s)-z(0) + -(z( S y-z(0y))}, s = q\ (A3) 



^(■s) = ^ — , r = v^f (v^+ - V T + ~ T -) > T± = (m B ±m K ) 



and the input given in Table IV The values fi(q = 0) stem from LCSR calculations. 
The form factors are shown in Fig. [9j extrapolated to high q 2 . Here, the agreement with 
preliminary lattice results by Liu et al. [26], which are shown as well, is good. 
Form factors stemming from different methods and parameterizations relevant to B — >■ 



K£ + £ decays at low recoil are compared in Fig. 10 Here we show / + and the ratio fo/f- 



entering B — >■ Kt + t decays divided by the corresponding (extrapolated) ones used in our 
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FIG. 10. Ratios of /+ and fo/ f+ from extrapolated LCSR by [37] (solid lines), extrapolated LCSR 
with simplified series expansion (38] (dashed), a relativistic quark model [39] (dash-dotted) and 
unquenched lattice calculations [26] (points), over the corresponding extrapolated form factors [7J 
used in this work. 

analysis, from Ref. [7]. Owing to a different g 2 -shape the extrapolated LCSR ones from [37] 
(solid lines) grow much larger towards very low recoil, and differ most strongly from the 
findings of [7J and lattice [26]. Both findings for / + in the relativistic quark model [39J and 
LCSR from [37J combined with a simplified series expansion [38] exhibit a shape similar to 
the ones from [7J . The differences between / + and f / f + from [7J and [3H1 EH] are within 



25% and 21%, respectively, which is within the uncertainties covered by Table IV 

The validity of the lowest order improved Isgur-Wise form factor relation Eq. (2.4) can be 

quantified by looking at deviations of 

q 2 fAq 2 ) 



R T (q 2 



m B (m B + m K ) f+(q 2 ) 
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FIG. 11. The relation Rx(q 2 ) given in Eq. (A4) from extrapolations of the LCSR form factors 
of [7J (blue band) and lattice [26] (red points). The grey band indicates the A/m^-correction to 



Eq. (2.4), see text. 



Form Factor f. t 


Resonance 


Si(0) 


b\ 


u 

So 
St 


m res ,+ = 5.412 GeV 

no pole 
m resT = 5.412 GeV 


n 34+0.05 
n 34+0.05 

u -°^-0.02 

n 30+0.05 


9 1+0.9 
z - i -1.6 

_4 3+0-8 

9 9+I.O 
z - z -2.0 



TABLE IV. Input to the B — >■ K form factor parameterization Eq. (A3) from [7J. 



from k, given in Eq. (2.5). As can be inferred from Fig. 11, for both the LCSR extrapola- 
tion [7j and lattice results [26J the agreement is good near the kinematical endpoint. For 
smaller dilepton masses is smaller than k 1. The agreement improves somewhat if 



the kinematical prefactor on the right-hand side of Eq. (A4) is replaced by one. Note that 
Rt{q 2 ) obtained with the form factors from 



behaves very similar to the one of [7J. 
Including a positive- valued l/m& HQET form factor 5+\q 2 ), as, for instance, from [TT] . 



the ratio St/S+ calculated from Eq. (2.3) would increase further. The form factors 8±^ are 



defined in terms of the heavy 6-quark field h v as [11J 

(K(k)\ sWpK \B(p)) = 5 ( °\q 2 )(p + k), + 5^\q 2 )q, . (A5) 

Lattice studies for B — > it [40J suggest that they are not larger than the estimate by dimen- 
sional analysis. 

We recall that the impact of St on the B — > K£ + £~ decay observables at low recoil is 
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subleading. This hampers on one side the extraction of fx from data, on the other side 
reduces the theoretical uncertainties. Ultimately it is desirable to know all form factors 
more precisely from lattice QCD, or other means. 



[1] C. Bobeth, G. Hiller and D. van Dyk, arXiv: 1105. 2659 [hep-ph], to appear in the Proceedings 

of DISCRETE2010 (La Sapienza, Rome), December 6th-llth 2010. 
[2] A. Ali, P. Ball, L. T. Handoko and G. Hiller, Phys. Rev. D 61, 074024 (2000) [arXiv:hep- 

ph/9910221]. 

[3] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 73, 092001 (2006) [arXiv:hep- 
ex/0604007]. 

[4] C. Bobeth, G. Hiller and G. Piranishvili, JHEP 0712, 040 (2007) [arXiv:0709.4174 [hep-ph]]. 
[5] M. Beneke, T. Feldmann and D. Seidel, Nucl. Phys. B 612 (2001) 25 [arXiv:hep-ph/0 106067]. 
[6] M. Beneke, T. Feldmann and D. Seidel, Eur. Phys. J. C 41 (2005) 173 [arXiv:hep-ph/04 12400]. 
[7] A. Khodjamirian, T. Mannel, A. A. Pivovarov and Y. M. Wang, JHEP 1009, 089 (2010) 

[arXiv: 1006.4945 [hep-ph]]. 
[8] B. Grinstein and D. Pirjol, Phys. Rev. D 70 (2004) 114005 [arXiv:hep-ph/0404250]. 
[9] M. Beylich, G. Buchalla and T. Feldmann, Eur. Phys. J. C71 (2011) 1635 [arXiv:1101.5118 

[hep-ph]]. 

[10] G. Buchalla and G. Isidori, Nucl. Phys. B525, 333-349 (1998) [hep-ph/9801456]. 

[11] B. Grinstein and D. Pirjol, Phys. Lett. B 533 (2002) 8 [arXiv:hep-ph/0201298]. 

[12] C. Bobeth, G. Hiller and D. van Dyk, JHEP 1007, 098 (2010) [arXiv: 1006.5013 [hep-ph]]. 

[13] C. Bobeth, G. Hiller and D. van Dyk, JHEP 1107, 067 (2011) [arXiv: 1105.0376 [hep-ph]]. 

[14] M. Patel, T. Blake for the [LHCb Collaboration], Conference report prepared for the 2011 

Europhysics Conference on High- Energy Physics, LHCb-CONF-201 1-038. 
[15] T. Blake, [arXiv:1109.6582 [hep-ex]]. 

[16] A. Golutvin, Talk given at La Thuile 2011, LHCb-TALK-2011-029. 
[17] T. Hurth and D. Wyler in J. L. Hewett et al, arXiv:hep-ph/0503261. 
[18] N. Isgur and M. B. Wise, Phys. Rev. D42, 2388-2391 (1990). 

[19] D. A. Demir, K. A. Olive, M. B. Voloshin, Phys. Rev. D66 (2002) 034015 [hep-ph/0204119]. 
[20] A. K. Alok, A. Dighe and S. Uma Sankar, Phys. Rev. D78, 114025 (2008) [arXiv:0810.3779 

21 



[hep-ph]]. 

[21] P. Z. Skands, JHEP 0101, 008 (2001) [arXiv:hep-ph/0010115]. 



[22] EOS Collaboration, a code for flavor observables http : //project .net .physik. 
tu-dortmund . de/ eos/ . 



[23] J. Charles et al. [ CKMfitter Group Collaboration ], Eur. Phys. J. C41, 1-131 (2005) [hep- 
ph/0406184]. We use the numerical results as presented at ICHEP10. 

[24] K. Nakamura et al. [Particle Data Group], J. Phys. G 37, 075021 (2010). 

[25] [Tevatron Electroweak Working Group and CDF Collaboration and DO Collab], 
arXiv:0903.2503 [hep-ex]. 

[26] Z. Liu, S. Meinel, A. Hart, R. R. Horgan, E. H. Muller and M. Wingate, [arXiv:1101.2726 
[hep-ph]]. 

[27] J. T. Wei et al. [BELLE Collaboration], Phys. Rev. Lett. 103, 171801 (2009) [arXiv:0904.0770 
[hep-ex]]. 

[28] T. Aaltonen et al. [CDF Collaboration] [arXiv:1107.3753 [hep-ex]]. 

[29] U. Egede, T. Hurth, J. Matias, M. Ramon and W. Reece, JHEP 0811, 032 (2008) 

[arXiv:0807.2589 [hep-ph]]. 
[30] J. Walsh, Talk given at the 46th Rencontres de Moriond: Electroweak Interactions and Unified 

Theories 2011, (La Thuile), March 13th-20th 2011. 
[31] T. Aaltonen et al. [CDF Collaboration], arXiv:1108.0695 [hep-ex]. 
[32] C. Wacker, Diploma Thesis, TU Dortmund, November 2011. 

[33] G. Hiller, arXiv:1106.1547 [hep-ph], to appear in the proceedings of the 46th Rencontres de 
Moriond: Electroweak Interactions and Unified Theories 2011, (La Thuile), March 13th-20th 
2011. 

[34] A. Ali, G. Kramer and G. -h. Zhu, Eur. Phys. J. C47, 625-641 (2006) [hep-ph/0601034]. 

[35] P. Ball and R. Zwicky, Phys. Rev. D71, 014029 (2005) [hep-ph/04 12079]. 

[36] W. Altmannshofer, P. Paradisi and D. Straub, arXiv:1111.1257 [hep-ph]. 

[37] P. Ball and R. Zwicky, Phys. Rev. D71, 014015 (2005) [hep-ph/0406232]. 

[38] A. Bharucha, T. Feldmann and M. Wick, JHEP 1009, 090 (2010) [arXiv:1004.3249 [hep-ph]]. 

[39] A. Faessler, T. Gutsche, M. A. Ivanov, J. G. Korner and V. E. Lyubovitskij , Eur. Phys. J. 

direct C4, 18 (2002) [hep-ph/0205287]. 
[40] E. Dalgic et al, Phys. Rev. D73, 074502 (2006) [hep-lat/0601021]. 



22 



